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We apply the recursive stochastic state selection method, vifhich is a new method for 
Monte Carlo study we have recently developed, to quantum spin systems with positive 
definite Hamiltonians. Through numerical studies of two-dimensional Ji — J2 Heisenberg 
model on a square lattice with unfrustrated couplings Ji = 1 and J2 = — 1 and with 
non-frustrated ones Ji — 1 and J2 = 0, we find that a kind of equilibrium is realized in 
these systems. We also observe that in this equilibrium we can obtain a quite accurate 
estimate of the energy eigenvalue for the system's ground state. Statistical relative 
errors in our results are 0.03% for the 36-site unfrustrated model and 0.06% for the 
64-site non-frustrated model. 
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1. Introduction 

Quite recently we have proposed a new method for Monte Carlo calculations to evaluate energy 
eigenvalues of quantum spin systems, which we call the stochastic state selection (SSS) method. 
This method enables us to obtain expectation values of powers of the Hamiltonian even when limited 
computer memory resources are available. In the SSS method we numerically select "active" basis 
states using random choice matrices. A random choice matrix is a diagonal matrix whose elements 
are random variables to follow two- valued probability functions, named on- off probability functions, 
which are defined by an initial trial state of the system. Applying the method we succeeded in 
obtaining reliable values for energy of the first excited states in the 64-site Shastry-Sutherland 
model. Further, we have modified the method to develop the recursive stochastic state selection 
(RSSS) method, where the on-off probability functions depend on the intermediate states of the 
system. '^^ A merit of the RSSS method is that errors of observed expectation values increase less 
rapidly than those in the SSS method when the power of the Hamiltonian rises. We applied the 
RSSS method to the anti-ferromagnetic Heisenberg spin one-half system on the triangular lattice, 
which is a typical strongly frustrated system to which other methods such as the ordinary Monte 
Carlo methods and the perturbative calculations are hardly applicable. Our result for the lowest 
energy eigenvalue of the 36-site system is within one percent of the exact value. 

In this paper we apply the RSSS method to quantum spin systems whose Hamiltonians are 
positive definite, namely whose Hamiltonian matrix elements are all non-negative. Tangible exam- 
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pies we employ here are the two-dimensional Ji — J2 Heisenberg model on a square lattice^' with 
unfrustrated or non-frustrated couplings. Numerical studies suggest that in these systems we are 
able to obtain, through repeated operations of the Hamiltonian and the random choice matrix, a 
normalized state which contains a finite and definite portion of the ground state. This implies a 
kind of equilibrium, which we call the RSSS equilibrium. Then we show that, by virtue of this equi- 
librium, one can easily estimate the energy eigenvalue of the ground state from the normalization 
factor. 

In the next section we briefly explain the RSSS method and give a definition of the RSSS 
equilibrium. Section 3 adds an extended analytical discussion on models with positive definite 
Hamiltonians, which indicate that they provide probable candidates to realize the RSSS equilibrium. 
In sections 4 and 5 we investigate the Ji — J2 model with Ji = 1 and J2 < 0. In section 4 we 
present numerical results to examine the discussions in sections 2 and 3. Here we demonstrate 
the realization of the RSSS equilibrium using small systems for which the exact ground states are 
calculable. First we study the non-frustrated model on a 4 x 4 square lattice setting J2 = — 1. Then 
the unfrustrated case ( J2 = 0) is investigated on the 6x4 lattice. Section 5 shows results on larger 
lattices, the 6x6 lattice with J2 = — 1 and the 8x8 one with J2 = 0. Here we describe how we 
obtain an estimate of the energy eigenvalue in the RSSS equilibrium. Our results thus obtained 
are —40.644 it 0.013 for the 36-site unfrustrated model, which is comparable to the value —40.659 
obtained from ref. 4, and — 43.099 ± 0.025 for the 64-site non-frustrated model, which is only 0.02% 
higher than — 43.107,^^ the most accurate value to our knowledge. The final section is devoted to 
summary and discussions. 

2. Recursive Stochastic State Selection Equilibrium 

First we briefly review the RSSS method."^) Let us denote the Hamiltonian of a system by H, a 
basis by {| i)}, the number of the basis states by A'^v and a trial function by | We introduce 

the random choice matrices M^j^(m)j = diag.jr/j^^^ 772"*^ ' ' ' ^VNy} (™ ~ 1)2, ••• ,L) in order to 
calculate 

E^^}{L) = I HM^CL)yHM^^^,,-^^ ■ ■ ■ HM^^^} I ^^°^), (1) 

where the random variable r]^"^^ is generated following the on-off probability function P^"^\r]). In 
the RSSS method we define 

using a positive parameter e and the coefficient cf^ in the normalized intermediate state 

I = ^ I i)c|™'~^^ which is proportional to HM^^(m-i)j ■ ■ ■ HM^^(i)y | Starting from a 
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given I ^(°)) we can recursively calculate 

where C^"^^ (> 0) is the normalization factor calculated by 
We also define | ^ira-i)-^g(m-i) 



(3) 



(4) 



(5) 



with the normalization condition (x^"*~^^ | x^^~^^) = 1- It should be kept in mind that for any 
state \ ^) = J2 \ with 6j's which are irrelevant to {Vi^^}, the statistical average of the inner 
product between | $) and | is zero, 



(6) 



because 



Jo 



(7) 



Note that, as was discussed in ref. 3, this average should be calculated with a fixed 
|^(m-2)|^ . . . jj^ ^j^jg papgi- -vvg call such states as | T^(™~^))(yr(™'~^) random 

states. For any random state | Q) which appears here we assume that the deviation of (<& | is 
negligibly small so that ($ | r^) ~ holds with no statistical average. Also remember that 



Jm-l) 



' 0<\4"'-'^\<e 



(m— 1) 



im-1) I 



E 



c. 



(m-l). 



where we use (7) and 



\c^-'^\<e 



(m-l) 



(8) 



(m) 



(9) 



Now we discuss on the RSSS equilibrium. Let | ^e) denote the exact eigenstate with the largest 
eigenvalue E. We divide the intermediate state | V^^"*)) into a part which is proportional to | ■0e) 
and the rest, 



(10) 
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where 



w 



(m) ^ 



= (V'E I V'^'")) , (11) 

and we request (C^™^ | C^™^) = 1- Note that 

(V'E I C^"'^)s^'"^ = , (13) 

and 



2 



+ 



2 



1 , (14) 



by definition. What we mean by the RSSS equihbrium is that there exists a hmit w^^"^^ defined by 

^ ms+mt — l 

hm — y = w^^i^ , (0 < w^^"^ < 1) , (15) 

m=ms 

where w^^'^^ is independent of mg, whenever mg is greater than or equal to some value of m, say 
mo- 

Let us derive a relation from this RSSS equilibrium for later use. Since from (3), (5) and (10) 

I ^(m+l)^ (jim+1) ^ HM^^^ra+i)y I V'^"'^) = ^ { I V'^"'^)+ I } 

= H { \ ^e)w'^"'^+ I C^""^) 5^""^+ I X^"'^)9^"'^ } 

= E I ^e)w^"'^ + H I C^"*^) s^"*^ + H I X^"'^)5^"'^ , (16) 

we obtain, using (11) and (13), 

^{m+l)fj{m+l) ^ I ^(m+l)^ ^ ^^{m) ^ ^^^^ | ^i.rn)^g{m) ^ (^7) 

which leads 

(m+l) ^ y;M (18) 

(^(m+1) ' ^ ' 

if the second term in the right-hand side of (17) is negligible. Then we can expect 

E = C(™+^) , (19) 

for sufficiently large values of m. Thus we become aware that the value of E can be estimated by 
the normalization factor C^'^'^^\ Actually, contributions from the second term in (17) should be 
taken into account when we consider the fiuctuation of C^'^'^^\ 
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3. RSSS Equilibrium for Positive Definite Hamiltonians 

In this section we deal with several quantities concerning to the systems with positive definite 
Hamiltonians, which will be referred to in our numerical work. After the manner of the previous 
section we make our discussion as simple as possible, neglecting all fluctuations in these quantities. 

In order to make our analysis concrete, we here limit ourselves to a case where 



V^(o)) =1 ^e) = J2\ i)fi ifi > for all i) . 



(20) 



We suppose we have chosen an adequate basis {| i)} so that all /j's as well as hij = {i \ H \ j)'s are 
non-negative, which is always possible for positive definite Hamiltonians. It should be noted that 
all coefficients in the expansion of | i^^^^), \ i'^'^^), ■ ■ ■ are also non- negative in this case, because in 



the relation between c|"^'* and c^"* ^\ 



(21) 



which we learn from (3), hij > for all i and j and rjj /C^'^^ > for all j by definition 
We pay our attention to a relation led from (16) which defines | ^("*+^))C^"*+^\ 



(j(m+l) 



w 



(in) 



W 



(m) 



+ I I 

+ 2£;2i(;(™)(V'E I X^'"^)5^"'^ + 2s('")(C('"^ | \ X^"'^)9^"'^ , (22) 

where we used (13) and (14). For sufficiently large value of m we expect that both (C^"*^ | \ C^"*^) 
and (x^"*^ I \ x^^^) independent of m, 



2C > 



2X 



(23) 
(24) 



where i?2C ^^'^ denote positive constants determined by e and H. We also expect both 
cross terms in (22) are negligible because statistical averages of them vanish.''^) Validity of these 
assumptions will be numerically examined in the next section. We thus obtain a relation 



^(m+l) 



w 



(m) 



+ H., 



2C 



w 



(ni) 



Let us examine then. Using (8) we obtain 



(m) 



(m) 



\Ci l<e 



14 '|<e 



(m) I 



5 ^E 



2X 



(m) I 



(m) 



(25) 



(26) 



The last equality follows from the fact that all c-™^ are non-negative here. Note that 



^\c^'>\<e 



< 1 because of the normalization condition. For m = 1 we obtain from (21) 



(27) 
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using the relation ^ hijfj = Efi which results from H \ i/j-e) = E \ ■0e)- Therefore 

i i i j 



(28) 



where we use (18), notifying that = 1 for (20). We assume that the second term in the right- 
hand side of (28) is negligible because sums over i and j in it would be enough to promote the 
cancellation in (t]^'^ — 1 ) keeping C*-^^ almost irrelevant to the sampling. Therefore we acquire a 



(29) 



relation e c^^^ ~ w^^^e /j. In the same manner we obtain 



,(1) 



which leads, taking (( {{(vf^rjl^^ - l)))fixed{)7(i)})) = ^ ^^^^ account, e^-cf^ ~ w^'^^^Y^ifi- Thus 
we expect 



(m) 



'-)) < w^-hY^n- y: [ 



(m) 



(30) 



\c\ '\<e 



holds for any m. In order to propose a relation to be expected in the RSSS equilibrium we add, 
based on the above discussions, the following assumption for sufficiently large m, 



|c- '\<e 



(m) 



~ K 



(m) 



(31) 



(32) 



where K (0 < K < 1) is a constant defined by e. We will see that these assumptions are acceptable 
for systems we numerically study. 

Now we reach to a function which determines C^"*"^^) from w^"^^ , 



Q{m+l) 



w 



(m) 



H2C I 1 



w 



(m) 



(Gw^"''> - 



w 



(m) 



= [E"" - H2c) 

Together with (18) this leads to a relation between w^"^~^^^ and w^'^\ 

Ew'^^^ 



w 



(m+l) 



^J {E^ - H2c) + GH2^w^-^) + {H2C - KH2^) 

Next notify it is mathematically guaranteed that when we define {xn} by 

Xn+l = f{Xn) , < Xo < 1 , 

with 

px 



(33) 



(34) 



V ax"^ + 26a; + c 



(0 < a; < 1) , 



(35) 



(36) 
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the value of Xn converges to a finite value between and 1 provided a>0,6>0,c>0,p>0 and 



p/\/a + 2b + c < 1. The reason for it is that, under these conditions, /(O) = 0, /(I) < 1 and 



f"{x) = -p 



2abx^ + (6^ + 3ac).T + 26c 



V ax"^ + 26a; + c 



< . 



(37) 



Let a = E"^ — H2(^, b = Gi?2x/2, c = i72C — -f^-f^2x and p = E then in (34). Since E is the largest 
eigenvalue of the positive definite H, it is clear that a > 0, 6 > and p > by definition. In addition, 
we see that a + 26 + c = {E^ - i^sc) + G'^2x + {H21; - KH2^) = E^ + iJac {G - K) > E^ = 
because we obtain G — K > Gw^"^^ — K ^ [5^'"-'] > using (32). Therefore (34) suggests that the 
RSSS equilibrium will be realized when i72C > KH2X, which seems mostly fulfilled irrespective of 
values of e.^^ 

Finally let us point out that we can obtain an equation for w^^'^^ combining (19) and (33), 
which is 

' ^ + Gi?2xw^'^^ - {E^ - H2C + KH2y)) = . (38) 



[E^ - i?2c) 
The relevant solution is 



w 



W 



(^1) = -q + ^/^TT+ 



K 



1 

2 ■ 'E^ 



GH. 



2X 



2X 



2C 



E^-H' 



2C 



(39) 



which should be compared with w^'^^ to examine the validity of (33). 



4. Numerical Study on Small Systems 

Here we numerically study two quantum spin one-half systems on small lattices in order to 
examine assumptions and relations discussed in previous sections. 

First we investigate the 4 x 4 Ji — J2 model with couplings Ji = 1 and J2 = —1. This system 
size is too small to satisfactorily suppress fluctuations of the observed data for e > 0.05, yet it has 
a merit that we can easily obtain the exact ground state by means of the Lanczos method. The 
Hamiltonian of the system is 

= ^Y^'^i-'^3 + ^^ '^i- '^3' ' (40) 
(nn) (nnn) 

where the first summation is over the nearest neighbor pairs and the second one over the next- 
nearest neighbor pairs on the square lattice. In the following calculations we employ, instead of 

H = li + U^ {-Hj,j,) U , (41) 

where / is the identity operator and U denotes a unitary transform to make off-diagonal elements 
hij {i 7^ j) non-negative. The value of the parameter I is determined so as to it ensures that all 
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diagonal elements of H are non-negative. We put I = 8 here. Let | ^e) be the eigenstate for the 
largest eigenvalue of H, which is expressed by as we did in the previous section. The number 
of basis states for | -^e) is 12, 870 and we obtain the exact value of E for this system is 26.633, 
which is in agreement with the result reported in ref 4. We denote the exact value of E by Eq 
hereafter. Note that, with this positive definite H, all coefficients in the expansion of | V'e) are 
non- negative. Starting | ■^e) we calculate | V^^"*)) up to m = 1000 using the /c-th random number 
sequence (1 < A; < nsmpi) to generate {q} = {{r/^'")}, {r/^"""^)}, • • • , {??^^^}}, which we refer to as a 
rns-sample hereafter. Although the initial transient process terminates up to m ~ 5, we need to 
carry out such long runs because some quantities measured with e > 0.05 show much fluctuations 
ranging a few hundreds of m. Figures 1-7 show results from one rns-sample. The number of basis 
states whose coefficients are non-zero is saturated to be ~ 5 x 10^ (~ 1.3 x 10^, ~ 3.8 x 10^) in 
the expansion of M^^(m)y \ V-^"""^^) for e = 0.1 (0.05, 0.02), while it is ~ 8.5 x 10^ (~ 1.2 x 10^, 
~ 12,870) in the expansion of HM^^^m)^ \ ip^'^"^^). The results are not much different for other 
rns-samples. 

Figures 1-4 are to numerically examine the assumptions employed when we obtain (25) from 

(22) . Figures 1, 2 and 3 plot quantities which appear in the right-hand side of (22) obtained 
for e = 0.1, 0.05 and 0.02, respectively. We see that cross terms 2E'^w^'^^ {ijfE \ x^^^)9^^^ and 
25(m)^^(m) I ^2 I ^im))gim) distribute around and are negligible compared to other three terms 
E^ [luM]^, (^M I I (M) (l- [u;M]^^ and (x^""^ | \ X^""^) [f^""^]^- Figure 4 shows 
((M I ^2 I ^(m)^ ^^(m) I ^2 | ^(m)^ calculated from | C^""^) defined in (12) and | x^""^) defined 
in (5). Although there are much fluctuations in data for (C^"*^ | \ C^^^), our assumption in 

(23) that they are nearly equals to a constant H21; seems acceptable at least when e = 0.02. As for 
^•^M I ^2 I T^M)^ we see the data are stable to support that (24), namely (x^"*^ I H'^ \ X^^^) — -f^2x5 
holds. Note that the condition > KH2y^ (0 < if < 1) is guaranteed from these data. 



In fig. 5 we plot ^| (m), cj™-* together with K which is estimated from the average of 
them with m = 800, ■■ ■ ,999. Figure 6 compares [5^"*^] calculated from (5) with Gw^'^^ — K 
{G = eY^fi), where we calculate w^'^^ from (11) and use K obtained in Fig. 5. These results 
suggest that the assumptions (31) and (32) are probably good for e = 0.05 and e = 0.02, although 
data for e = 0.05 fluctuates a little and Gw^"^^ — K for e = 0.02 slightly overestimates [y*^™-*]^. 
Results for e = 0.1, on the contrary, are not fully convincing because of their fluctuations. In 
Fig. 7 we present data for w^'^^ together with the solution w^^'^^ in (39), where we substitute the 

r "1 2 

average of ^, cj™^ , (C^"*) | \ C^""^) and (x^"*) | \ X^""^) over m = 800, • • • ,999 for 

K, H2(^ and i?2xj respectively. We see the agreement is satisfying. In order to make a numerical 
examination of (15) we calculate, changing mg, averages of w^'^^ 's for finite values of mt- The results 
prove to be almost independent of mg. Although data with e = 0.1 give us a little fluctuating 
results such as 0.575 (mg = 200, mt = 800) and 0.584 (mg = 800, mt = 200), others are very 
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Stable. They are, for instance, 0.7885 (mg = 200, mt = 800), 0.7883 (mg = 600, m-t = 400), 
0.7885 (ms = 600, mt = 200) and 0.7881 (mg = 800, mt = 200) with e = 0.05 and 0.9415 
{rris = 200, mt = 800), 0.9416 (mg = 600, mt = 400), 0.9417 (mg = 600, mt = 200) and 0.9415 
(mg = 800, mt = 200) with e = 0.02. Prom these results we conclude that the RSSS equilibrium is 
established in this system for all values of e, and that our discussion to lead (39) is valid. 
Figure 8 plots 

r'-'K,,- — Ec«. (42) 

where C^^y^ denotes the normalization factor calculated from (4) for the A;-th rns-sample. The 
total number of these samples is rigmpi = 100. The error is estimated by 



= 2 



^smpl 1 

The exact value of E is also shown in the Figure by a dashed line. We observe that data for e = 0.02 
support the approximate relation 

E = ((C("^+^))) , (44) 

which follows from (19). Data with e = 0.05 and e = 0.1 fluctuate much more than those with 
e = 0.02, but we can still see they are close to the exact value of E. These fluctuations of 
((C^"*"*"^^)) would be mainly ascribed to the fact that the second term in (17), which we evaluate to 
be KV-E I x^"'^)9^'"^/w^"'M < 0.03 (0.07, 0.14) for e = 0.02 (0.05, 0.1) in our measurement, is not 
completely negligible. We think this term is responsible for the e-dependency of ((C*^™"*"^))) which 
will be discussed in the following section. 

Now let us show results for the J2 = case studied on a 6 x 4 square lattice. Here we employ 
Z = 6 in (41), for which £'0, the exact value of E, becomes 22.553. Although it is more time- 
consuming to investigate this lattice size, we can obtain results which fluctuate less compared to 
those for the 16-site J2 = — 1 model. The number of the basis states in the the exact eigenstate is 
2, 704, 157. The number of basis states whose coefficients are non-zero is ~ 2.4 x 10^ (~ 4.5 x 10^) 
in the expansion of M^^(m)y | V'^™"^^) for e = 0.05 (0.01), while it is ~ 6.5 x 10^ (~ 8.6 x 10^) in 
the expansion of HM^jj(m.)y \ 

Figure 9 presents data of [yMj^ with e = 0.05 and e = 0.01 for one rns-sample, together with 
Gw^'^^ — 1, where we substitute 1 for ii' . For e = 0.05 it is difficult, due to fluctuations in w^"^\ to 
see that [5^'"^]' and Gw^"^'' — 1 are comparable. When we take the less value e = 0.01, however, we 
can ensure that ~ GiuM - 1. In Fig. 10 we plot lyM. We see it is quite low (< 0.2) when 

e = 0.05, while it increases to ~ 0.6 with e = 0.01. Figure 11 shows ((CM)) over 100 rns-samples, 
which should be compared with Eq plotted by a dashed line. We see that when e = 0.01 values of 



10 Tomo MUNEHISA and Yasuko Munehisa 

((C^"*))) after the transient decrease are nearly equal to Eq. In addition, it is noticeable that data 
with e = 0.05, whose overlaps with the exact eigenstate {w^"^^) are less than 0.2, are also closely 
located to Eo when m > 10. 



5. Numerical Results on Large Systems 



In this section we study the Ji — J2 model on large lattices, the 36-sitc J2 = — 1 model with 
i = 18 in (41) and the 64-site J2 = model with / = 0. Here we measure only C^"^^ and 



m-l) 



M 



{r;(™)} 



(m-l) 



(45) 



for m = 1, 2, • • • , mmax- The quantity (45) is measured to make sure that [5^"^"^^] does not diverge 

for large vahics of m. 

First we show results for the 36-site J2 = — 1 model. Since we do not know the exact eigenstate 
of this system, we employ the Nccl state as the initial trial state | ip^^^). For a few hundred of 
rns-samples we calculate C^"*) and [c/^"*"^)] with 

n^ma.x — 1000 with several values of e between 
0.02 and 0.1. For any of the rns-sample we observe [5'^'"-*]^ '-^ 6.3 when m > mo = 200 regardless 
of the value of e. The upper bound of the number of basis states whose coefficients are non- 
zero ranges between ~ 8 x 10^ (for e = 0.1) and ~ 1.9 x 10^ (for e = 0.02) in the expansion of 



M 



m— 1) \ 



, while it is within ~ 6 x lO'' (for e = 0.1) and ~ 1.4 x 10^ (for e = 0.02) in the 



expansion of i^M|^^(m)} | V*-™^^"*)- Figure 12 and 13 plot {{C^"^^))n^^pi with rismpi = 100 for several 
values of e. Dashed lines in these Figures present the value of E obtained from ref. 4, which amounts 
Eq = 58.659 for this system. We see that the data for m > uiq are near Eq with some fluctuations. 
Let us denote the normalization factor for | ^(™)) obtained in the A;-th sample by C^^^ as before. 
In order to estimate the exact eigenvalue of E based on (19) with a better precision, we take an 
average of C'j^'*^ not only over k but also over m = mo, mo + Am, mo + 2 Am, • • • , with some 
positive integer Am. Namely, we calculate 

"smpl rtmax 



1 



(ni, 

for 7 = 1, 2 and 



+ 1) 



2^ 2^ rw. 



k=l n=0 



m-max - "'0 

Am 



(46) 



Er = 2\ 



C2 



(47) 



(^max "I" 1) ■ ''^smpl 1 

It should be noted, however, that if Am is not large enough there will be a correlation between 
C^y^ and in contrast to that C^y^ and C^y \ are statistically independent by definition 

iik ^ k' . We therefore determine Am by the two-sided 5% t-test of the hypothesis that there is no 
correlation between C'j^''^ (m > mo) and C'j^^^'"^- The result of the test suggests that Am should 
be more than 90, so we set Am = 100 for this model. The results thus obtained with ngmpi = 100 
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(200) for e < 0.07 (> 0.07) are shown in Fig. 14. We observe that they are located close to Eo but 
small discrepancies, which seem to be linearly dependent of e, still remain. So we carry out the 
weighted fit by the method of lease squares, assuming that 

C = Co-Ae, (48) 

where Cq and A are positive constants to be pursued. Results of the fit using data with 0.02 < 
e < 0.1 (0.02 < e < 0.06) are Co = 58.644 ± 0.013 (Co = 58.673 ± 0.014). We plot the fitted line 
as well as the error of mean square for Co, which we denote by AC, by solid (dashed) lines in the 
Figure. These values lead, since / = 18 and Eq = 58.659 for this system, |ACo/(Co - 01 = 0.00032 
(0.00034), or \{{Eo - I) - (Co - 0}/(^o - 01 = 0.00037 (0.00034), in the fit using data between 
e = 0.02 and 0.1 (0.06). It is thus evident in this model that Co gives a good estimate for E, which 
strongly support the realization of the RSSS equilibrium in this large-sized system. 

Similar investigations are made on the I = 0, 64-site J2 = model, for which the most accurate 
value of E is known to be Eq = 43.107.^^ From several tens of rns-samples with mmax = 300 or 
JTimax = 500 starting from the approximate state used in ref. 1 we obtain the following results. 
The RSSS equilibrium seems to be realized for m > mo = 200, where [5^"*^]^ ~ 32 for all values 
of e we measured. The number of basis states whose coefficients are non-zero is, for m > nio, 
~ 1.3 X lO'' (~ 3.5 X 10^, ~ 7.9 x 10^, ~ 3.1 x 10^) in the expansion of M|^(™)} | tp^-""'^)) for e = 0.05 
(0.03, 0.02, 0.01), while it is ~ 9.5 x 10^ (~ 2.6 x 10^, ~ 5.9 x 10^, ~ 2.3 x 10^) in the expansion 
of HM^^(m)y I V'^"*"^^)- From the two-sided 5% t-test we conclude that Am = 20 is enough for 
^Wh^^^ to have no correlation with C^^^ (m > mo). Figure 15 plots values of C obtained from 
(f^max + 1) • f^smpi data, namely 6 x 100, 16 x 10, 6 x 50 and 6 x 36 data with e = 0.05, 0.03, 
0.02 and 0.01, respectively. Results from the least square fit to (48) are also shown in the Figure. 
The fitted value of Co is 43.099 it 0.025, which agrees well with Eq. With these values we obtain 
|ACo/Co| = 0.00058, or {{Eq - Co)/Eo\ = 0.00019. This provides another example to indicate that 
the RSSS equilibrium is established and that our way to estimate E in the equilibrium is quite 
powerful. 

6. Summary and Discussions 

In this paper we numerically study quantum spin systems with positive definite Hamiltonians 
by means of the RSSS method. We find that a kind of equilibrium, which wc call the RSSS 
equilibrium, exists. Wc also notify that in this equilibrium we can effectively estimate the energy 
eigenvalue of the ground state. 

In the RSSS method we recursively calculate the m-th normalized intermediate state from the 
(m — l)-th one, starting from an initial trial state. The procedure is as follows. First we operate the 
m-th random choice matrix to the (m — l)-th state, which drastically reduces the effective size of 
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the vector space. Then we successively operate the Hamiltonian, which again increases the number 
of the basis states relevant to the state. Finally we normalize the resultant state to obtain the 
m-th normalization factor and the m-th normalized intermediate state. What we assert is that, 
after repeating this procedure many times, the m-th state comes to contain a finite portion of the 
ground state which is irrelevant to m. This means that the system is in the RSSS equilibrium. 

Our results on the 16-site and 36-site Ji — J2 model with couplings Ji = 1 and J2 = — 1 as 
well as the 24-site and 64-site model with Ji = 1 and J2 = afford abundant evidence of the 
RSSS equilibrium in these systems. In addition, from the normalization factors stated above we 
obtain satisfying results which estimate the ground state energy on large lattices within 0.04% of 
precision. Our study on small and large lattices with various values of the parameter e suggests 
that the RSSS equilibrium is observable for any value of e, as far as it keeps the random choice 
matrices non-zero. 

A few remarks are in order. In the present study we are solely concerned with systems hav- 
ing positive definite Hamiltonians, to which many useful methods are known to calculate their 
energy eigenvalues. In fact, a large amount of numerical work for these systems has been already 
reported.^' What is the merit of our way to estimate the energy eigenvalue, then? A remarkable 
feature is that it is very simple, having little to do with details of the system such as the dimension- 
ality and so on. Even for more complicated Hamiltonians our way to study the RSSS equilibrium 
is as simple as for those studied here. 

How shall we think about the systems whose Hamiltonians are not positive definite? We 
suppose they also realize the RSSS equilibrium when the parameter e is small. Preliminary results 
on the frustrated Ji — J2 model and the triangular Heisenberg model are promising. 

There remain some fascinating tasks. The most important one is to give a rigorous proof of 
the RSSS equilibrium from mathematical point of view. Another task is to analytically evaluate 
fluctuations we neglected in this paper. For this purpose it is necessary to investigate e-dependency 
of each dropped term in (17) and (22). We expect that contributions from these terms result in the 
linear relation (48) for sufficiently small values of e. Arguments how to extract physical quantities 
other than the energy eigenvalues in the equilibrium would be also necessary. 
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Fig. 1. Terms in (22) for the J2 = —1 model with e = 0.1. Open circles are , filled diamonds are 

(x'""^ I I x'""^) filled squares are (C'"") | \ C^""^) open triangles are 2S2w("'){V'e | X^"'^)>S^"'' 

and filled triangles are 2s('"'{Cf'"' | \ x''"')^*'"'- 
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Fig. 2. Terms in (22) for the 16-site J2 = —1 model with e = 0.05, where the same symbols as those in Fig. 1 are 
used. 
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Fig. 3. Terms in (22) for the 16-site J2 = — 1 model with e = 0.02, where the same symbols as those in Fig. 1 are 
used. 
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Fig. 4. (C^™' I I C*™'} for the 16-sitc J2 = -1 model with e = 0.1 (open triangles), e = 0.05 (asterisks), e = 0.02 
(filled triangles) and (x*™' | H"^ \ X*™') with e = 0.1 (open circles), e = 0.05 (open squares), e = 0.02 (filled 
diamonds). Note that (C^"'^ | H"^ \ C'"') > (x*'"^ I \ x^"'^) always holds. 
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Fig. 5. X]|c("')|<e ["^i*"'] f*-"" 16-site J2 — —1 model with e = 0.1 (open circles), e = 0.05 (asterisks) and e = 0.02 
(open diamonds). We also show the value of K estimated from last 200 data with e = 0.1 (0.05, 0.02), which is 
K = 0.919 (0.849, 0.601), by a dotted (dashed, dot-dashed) line, respectively. 
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Fig. 6. [5'™']^ for the 16-site J2 = —1 model with e = 0.1 (open circles), e = 0.05 (open squares), e = 0.02 (open 
diamonds). We also plot Gw'"' - K with e = 0.1 (filled circles), e = 0.05 (filled squares) and e = 0.02 (filled 
diamonds), where G/e = X] /» = 82.08 and we use values of K estimated in Fig. 5. 
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Fig. 7. w''"' defined by (11) for the 16-site J2 = —1 model with e = 0.1 (open circles), e = 0.05 (filled squares) and 
e = 0.02 (open diamonds). Estimated solutions w''^''' for (39) (a dotted fine for e = 0.1, a dashed line for e = 0.05 
and a dot-dashed line for e = 0.02) are also shown for comparison. 
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Fig. 8. ((C'^'"'))nsmpi ("smpi = 100) for the 16-site J2 = —1 model with e = 0.1 (open circles), e = 0.05 (asterisks) 
and e = 0.02 (filled diamonds). Errors shown in the figure are Br^"*' defined by (43). The dotted line indicates 
Eo, the exact value of E. 
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Fig. 9. [5'"*']^ for the 24-sitc J2 — model with e = 0.05 (open circles), e — 0.01 (open diamonds). We also plot 
Gw<"'> - 1 with e = 0.05 (filled circles), e = 0.01 (filled diamonds), where G/e = 792.1, for comparison. 
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Fig. 10. w*™' defined by (11) for the the 24-site J2 = model with e = 0.05 (filled circles), e = 0.01 (open diamonds), 
together with estimated solutions w'*"^' for (39) represented by a dotted line (e = 0.05) and a dashed line (e = 0.01). 
Here we put if = 1 for both values of e. 
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Fig. 11. ({C''"'))„^„,p, (rismpi = 100) for the the 24-site J2 = model with e = 0.05 (asterisks) and e = 0.01 (filled 
diamonds). {(C'-^^})n,^pi = 24.01 when e = 0.05. Errors shown in the figure are Er^"'^ defined by (43). The dotted 
line indicates the exact value of E. 
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Fig. 12. ((C('"'))„,^pi (n^mpi = 200) for the 36-site J2 = -1 model for m < 100 with e = 0.1 (open circles), e = 0.05 

(asterisks) and e = 0.02 (filled diamonds). The initial trial state is the the Nccl state. Errors Er^"^^ defined by 
(43) are within symbols when e < 0.05. A dotted line indicates Eq, the value of E obtained from ref. 4. 
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Fig. 13. ((C'^'"'))na„pi (^^smpl = 100) for the 36-site J2 = — 1 model up to m = 1000 calculated under the same 
conditions as Fig. 12, with the same symbols as those in Fig. 12. We show Eq by a dotted line. 
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Fig. 14. C versus e calculated for the the 36-site J2 = — 1 model by (46) with mo = 200, mmax = 1000 and Am = 100. 

Errors are defined by (47). The number of Cj^j'^ 's used to calculate each C is 900 (1800) when e < 0.07 (> 0.07). 
The dotted line indicates Eq. The solid (dashed) line presents the result for the linear fit from the data with 
0.02 < e < 0.1 (0.02 < e < 0.06) obtained by the method of least squares assuming (48). Error bars at e = show 
the error of mean square for the fitted Co- 
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Fig. 15. C for the the 64-site J2 = model versus e. Wo set mo = 200, Am — 20 in (46). Total number of data 
to calculate C is 600 (160, 300, 216) when e = 0.05 (0.03,0.02, 0.01). Errors defined by (47) are also shown. The 
dotted line indicates the value of Eq. The solid line presents the result for the least square fit by (48) with the 
error bar to show the error of mean square for the fitted Co. 



